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PREFACE 

Ultraviolet  sensors  for  missile  defense  must  be  developed 
based  on  accurate  knowledge  of  the  radiative  environment  or 
background  within  which  such  sensors  will  be  used.  This  report 
represents  a  beginning  effort  to  provide  such  accurate  knowledge 
by  presenting  the  results  of  model  calculations  of  the  expected 
backgrounds  using  the  current  understanding  of  upper  atmospheric 
and  ionospheric  physics.  The  effort  is  a  part  of  the  Ultraviolet 
Phenomenology  work  supported  by  SDIO/T-SN,  Dr.  Barry  Katz, 
through  Dr.  Robert  E.  Huffman  at  AFGL. 

These  background  levels  are  suitable  for  use  in  planning 
ultraviolet  backgrounds  measurements  on  SDIO  observational 
programs  such  as  STARLAB,  IBSS,  etc.  Along  with  currently 
available  material,  this  information  has  been  made  available  to 
the  mentioned  SDIO  programs  through  Experimenter's  Working  Groups 
and  more  widely  through  WPD  32.1  Reviews  and  Phenomenology 
Steering  and  Analysis  Group  (PSAG)  meetings. 

Model  calculations  such  as  these  will  also  be  used  in  the 
critical  data  analysis  and  review  planned  for  new  observational 
data.  The  ultimate  goal  is  to  provide  the  best  available  models 
of  the  backgrounds  and  its  variability  to  SDIO  users  for  scene 
generation,  systems  architecture,  and  sensor  development 
applications. 

Until  new  data  are  available,  theoretical  treatments  of 
regions  with  complicated  spectral  structure  will  continue.  For 
this  purpose,  use  continues  to  be  made  of  existing  AFGL  data 
bases  from  the  AIRS  UV  imager  on  the  Polar  BEAR  satellite  and  the 
S3-4  VUV  Backgrounds  satellite  experiment  and  other  existing 
data . 
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1.  INTRODUCTION 


Computational  Physics,  Inc.,  working  jointly  with  the  Upper 
Atmospheric  Physics  Branch  of  the  Space  Science  Division  of  NRL, 
has  developed  optical  background  models  for  the  dayglow  and 
Rayleigh  scattering  of  sunlight.  The  wavelength  range  presently 
goes  from  about  800  A  to  10,000  A  (this  range  for  the  dayglow  and 
2010  to  8000  A  for  Rayleigh  scattering) .  Key  quantities 
calculated  by  the  models  are  synthetic  spectra  in  both 
megaRayleighs  A'‘  and  uWatts  cm"^  sr"’  n  '  for  any  observing  altitude 
and  look  direction  of  interest.  From  these,  we  further  obtain 
limb  profiles  by  integrating  these  spectra  over  specified 
wavelength  and  look  angle  intervals  and  then  joining  the  two 
components  (dayglow  and  Rayleigh  scattering)  together.  The 
integrations  allow  for  convolution  with  instrument  response 
functions  so  that  direct  comparisons  can  be  made  with 
observations.  This  report  will  describe  the  various  parts  of  the 
models  and  then  present  results  demonstrating  present 
capabilities.  We  refer  the  reader  in  advance  to  Table  1  which 
summarizes  observing  and  physical  conditions  being  addressed. 
(Section  5  is  dedicated  to  this  topic.) 

2.  DAYGLOW  MODEL 

Our  dayglow  model  provides  first  principles  calculations  of 
photoelectron  transport,  chemistry  and  photon  transport. 
Calculated  quantities  of  interest  are; 

1.  photoelectron  spectra 

2.  volume  excitation  rates 

3 .  ion  and  neutral  densities 

4.  intensities. 

Intensities  are  provided  by  codes  LIMB  and  LIMB1304.  LIMB  is  a 
major  code  developed  over  the  past  three  years  which  calculates 
intensities  for  a  variety  of  features.  The  geometry  is  spherical 
and  the  observing  point  may  be  either  above  or  within  the 
emission  region.  Intensities  are  specified  as  functions  of  look 
angle  which  may  range  from  0°  (zenith)  to  180°.  The  code 
contains  a  "smart"  integrator  which  examines  the  integrand  along 
the  line-of-sight ,  sets  up  integration  intervals  and  then 
performs  Gaussian  integrations.  Features  and  systems  currently 
handled  by  LIMB  include  the  O  lines  at  1356  A,  2972  A,  5577  A, 
6300  A,  7774  A  and  8446  A,  the  O*  834  A  triplet,  other  lines 
such  as  NI  1493,  O*  2471  A  and  O*  7320  A  and  finally  the  band 
systems  N^  LBH,  N^  VK,  N^  2PG,  N^  IPG,  NO  y  and  IN.  The 
volume  emission  rates  for  01  1356  A  and  0*  834  A  require 
special  attention  because  of  multiple  scattering.  The  needed 
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transport  calculations  are  performed  by  code  DSDA  prior  to 
running  LIMB. 

LIMB1304  is  also  a  major  code  recently  developed  to 
calculate  01  1304  A  limb  profiles  from  source  functions  obtained 
by  a  rigorous  frequency  redistribution  description.^*"  Profiles 
are  obtained  for  each  of  the  triplet  lines  and  for  both  sources 
(solar  1304  A  photons  and  photoelectron  impact  excitation) . 

3 .  RAYLEIGH  SCATTERING  MODEL 

Our  Rayleigh  scattering  model  provides  a  first-principles 
calculation  of  multiple  scattering  from  2010  A  to  8000  A.  The 
model  accounts  for  pure  absorption  by  O^  and  O3  and  allows 
for  a  range  of  surface  conditions  (different  albedos  and  surface 
heights) .  Rayleigh  scattered  sunlight  begins  to  appear  above  the 
dayglow  background  around  1750  A.  Below  this  wavelength,  pure 
absorption  by  Oj  prevents  the  sun's  radiation  from  penetrating 
deep  enough  into  the  atmosphere  to  produce  any  measurable 
backscatter.  Between  about  1750  and  2000  A,  the  backscattered 
spectrum  rises  about  three  orders  of  magnitude  and  has  a  high 
degree  of  structure  caused  by  O2  Schumann-Runge  band 
absorption.  Our  model  does  not  presently  address  this  region 
because  of  the  severe  computational  difficulties  associated  with 
the  structure.  Work  is  underway,  however,  to  provide  estimates 
of  radiance  limits  over  selected  band  intervals  within  this 
region. 

The  model  is  based  on  an  isotropic  formulation  that  performs 
a  one-time  matrix  inversion  to  determine  the  radiation  field  for 
all  angles  of  interest  to  the  external  photon  source  and  for  a 

given  model  atmosphere.*^'  For  Rayleigh  scattering,  the  results 
agree  to  high  accuracy  with  a  detailed  Monte  Carlo  model  of  the 
molecular  scattering.*^'  When  aerosols  are  present,  the  model 
reproduces  the  Monte  Carlo  results  to  10%  or  better.*^'  The  model 
assumes  plane  parallel  geometry  so  its  applicability  is  limited 
to  solar  zenith  angles  of  less  than  approximately  95°.  Tests 
indicate,  however,  that  up  to  95°,  the  results  are  accurate.*^’ 

To  examine  azimuthal  effects  referenced  to  the  sun 
direction,  the  radiation  field  is  separated  into  a  singly 
scattered  (SS)  and  multiply  scattered  (MS)  component.  The  MS 
component  is  treated  as  isotropic,  which  is  a  good  approximation 
for  monochromatic  scattering  and  the  SS  component  is  scaled  by 
the  Rayleigh  phase  function  appropriate  to  the  observing 
conditions.  In  the  presence  of  aerosols  or  clouds,  the  Mie  phase 
function  is  used.  The  present  report  addresses  Rayleigh 
scattering  only.  At  wavelengths  in  excess  of  320  nm,  the 
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effects  of  aerosols  and  clouds  can  be  important  and  may  be 
considered  in  a  future  report. 

Volume  emission  rates  calculated  by  the  above  model  are  read 
into  code  RAYLIMB  which  then  calculates  limb  radiances  in  units 
of  either  megaRayleighs  A'*  or  jiWatts  cm'^  sr'‘.  The  possible 
range  of  tangent  altitudes  is  from  approximately  100  km  to  -6370 
km  (-6370  refers  to  nadir  viewing) -  RAYLIMB  has  been  developed 
over  the  past  year  and  has  its  origins  in  code  LIMB. 

4.  SPECTRAL  MODELS 

The  outputs  from  both  of  the  above  described  models  are 
respectively  read  into  codes  SYNSPEC  and  RAYSPEC.  SYNSPEC  is  a 
major  code  written  in  IDL  which  generates  normalized  synthetic 
spectra  for  the  above  mentioned  band  systems,  scales  them  by  the 
calculated  total  system  limb  intensities  and  adds  all  such  scaled 
spectra  on  a  common  wavelength  grid.  For  LBH,  pure  absorption  by 
Oj  is  correctly  accounted  for  as  a  function  of  wavelength. 
Emission  lines  are  added  at  their  appropriate  wavelengths  after 
smoothing  which  is  consistent  with  the  chosen  resolution. 

Routines  of  R.  Conway  are  employed  to  obtain  the  noirmalized 
spectra  for  LBH,  2PG  and  VK.  These  routines  return  the  spectra 
at  the  user  specified  resolution,  rotational  temperature  and 
wavelength  grid  spacing.  A  routine  provided  by  R.  McCoy  is  used 
to  generate  the  normalized  spectrum  for  the  NO  y  bands.  The 
spectrum  contains  the  same  level  of  detail  as  those  of  Conway. 

The  remaining  systems  are  given  a  simpler  description.  Each  such 
spectrum  is  treated  as  a  series  of  lines  in  place  of  its  bands 
and  appropriately  smoothed  for  the  specified  resolution. 

RAYSPEC  is  also  written  in  IDL  but  is  a  much  simpler  code. 
The  wavelength  grid  is  composed  of  those  wavelengths  for  which 
limb  profiles  were  obtain  by  RAYLIMB.  It  basically  takes  the 
family  of  limb  profiles  calculated  by  RAYLIMB,  places  them  into  a 
two  dimensional  array  (versus  wavelength  and  tangent  altitude) 
and  proceeds  to  plot  spectra  for  user  specified  tangent 
altitudes. 

Important  improvements  have  recently  been  made  to  both 
SYNSPEC  and  RAYSPEC.  It  is  now  possible  to  integrate  the  spectra 
over  wavelength  and  look  angle  to  obtain  limb  profiles  for  either 
simulation  of  or  comparisons  with  measurements.  Many  wavelength 
intervals  are  permitted  within  each  run  for  which  the  user 
specifies  beginning  and  ending  wavelengths.  The  integration  over 
look  angle  gives  f ield-of-view  (FOV)  averaged  profiles  for  either 
a  circular  or  rectangular  FOV.  Finally,  the  resulting  profiles 
from  the  two  codes  may  be  added  together.  Results  will  follow 


showing  the  combined  profiles  for  several  integrated  wavelength 
intervals  and  three  circular  FOVs. 

5.  PHYSICAL  AND  OBSERVING  CONDITIONS 

The  dayglow  results  to  follow  are  sensitive  to  solar 
activity.  At  tangent  altitudes  such  as  300  to  400  km,  the 
dayglow  can  increase  by  more  than  an  order  of  magnitude  going 
from  low  to  high  activity.  This  is  due  to  the  increase  in  the 
photoelectron  population  which  is  the  main  source  of  the  dayglow 
and  to  increased  atmospheric  densities  accompanying  the 
increased  exospheric  temperature.  At  lower  tangent  altitudes, 
the  increase  will  be  on  the  order  of  two  to  three.  The  parameter 
which  specifies  solar  activity  in  our  model  is  the  10.7  cm  flux 
designated  Fig^.  Results  will  follow  for  values  of  70  and  250 
which  reflect  low  and  high  activity.  This  parameter  is  used  in 
an  empirical  relationship  applied  to  each  of  the  800  or  so  lines 
in  our  solar  spectrum.  The  relationship  and  the  base  spectrum 
come  from  Hinteregger.'^'  F,q7  is  also  a  parameter  appearing  the 

MSIS  83  model'"*  which  is  providing  our  model  atmospheres  in  this 
work.  [Two  parameters  are  actually  called  for:  an  averaged 
value  for  the  previous  three  months  and  a  value  for  the  previous 
day.  ] 


The  solar  spectrum  longward  of  2000  A  is  not  sensitive  to 
solar  activity  and  consequently  the  same  one  is  used  in  our 
Rayleigh  scattering  calculations  for  low  and  high  activity."^' 
Results  have  been  obtained  for  two  solar  zenith  angles  (SZA) : 

30’’  and  60°.  The  models  presently  handle  angles  up  to  about 
100°  for  photoelectron  produced  emissions  and  90°  for  Rayleigh 
scattering.  The  limit  for  Rayleigh  scattering  will  likely  be 
extended  in  future  work.  We  will  be  presenting  results  for  only 
the  first  of  the  two  SZAs.  This  is  because  the  background  does 
not  change  significantly  in  going  from  30°  to  50°  and  because 
of  the  large  number  of  figures  involved. 

Sixteen  adjacent  wavelength  intervals  have  been  selected  for 
generating  limb  profiles.  Each  is  100  A  wide  with  the  first 
starting  at  1400  A  and  the  last  ending  at  3000  A.  Tangent 
altitude  resolutions  of  1,  3  and  10  km  will  be  considered.  This 
parameter  refers  to  the  interval  seen  on  the  limb  and  is 
determined  by  the  satellite's  altitude  and  FOV. 

A  resolution  of  10  A  has  been  chosen  for  generating  dayglow 
spectra.  They  have  been  placed  on  a  uniform  grid  with  a  spacing 
of  2  A.  The  spectra  to  follow  begin  at  1200  A  and  include 
Hydrogen  Lyman  Alpha  at  1216  A  and  01  1304  A.  Their  limb 
profiles  were  calculated  specifically  for  inclusion  in  this 
report  but  should  be  viewed  in  a  somewhat  qualitative  sense.  In 
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actuality,  they  vary  with  model  atmosphere,  solar  activity, 
satellite  altitude  and  SZA.  The  given  results  apply  to  small 
SZAs  and  to  solar  activity  in  the  moderate  range. 

The  limb  profiles  in  this  report  will  be  given  versus 
tangent  altitude  (300  to  0  km)  rather  than  look  direction.  The 
advantage  of  this  is  that  the  profile  is  invariant  with  satellite 
altitude  provided  that  the  satellite  is  effectively  above  the 
emitting  region.  This  will  generally  be  the  case  except  for  HI 
1216  A  and  possibly  1304  A.  For  the  latter  feature,  the 
satellite  altitude  must  be  approximately  800  km  or  greater  to  be 
above  the  emitting  region. 

The  following  table  summarizes  the  parameters  discussed  in 
this  section. 


Table  1.  Parameters  related  to  physical  and  observing  conditions 


PARAMETER 

VALUE  or  REFERENCE 

Solar  spectrum 

References  9  &  10 

(X-rays  to  1000  A) 

Solar  spectrum 

Reference  12 

(2000  -  3000  A) 

Model  Atmosphere 

Reference  11 

F10.7  (10'^^  Watts  cm'^  Hz'‘) 

70  and  250 

Solar  zenith  angles 

30"  and  60° 

Wavelength  intervals 

Widths 

100  A 

Number 

16 

First 

1400  -  1500  A 

Last 

2900  -  3000  A 

Resolution  on  limb 

3,  10  and  30  km 

Spectral  resolution  (dayglow) 

10  A 

Grid  spacing  (dayglow) 

2  A 

Range  of  tangent  altitudes 

300  to  0  km 

6 .  VOLUME  EMISSION  RATES 

The  key  inputs  to  the  limb  codes  discussed  above  are  volume 
emission  rates.  In  the  case  of  the  dayglow,  the  rates  apply  to 
discrete  features,  these  being  either  atomic  lines  or  systems  of 
molecular  bands.  In  the  case  of  Rayleigh  scattering,  we  are 
dealing  effectively  with  continuum  emission  and  thus  may  define 
the  emission  rate  at  any  wavelength.  Selected  rates  will  follow 
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for  both  the  dayglow  and  Rayleigh  scattering  to  illustrate  the 
striking  contrast  in  altitude  dependence  which  can  occur  from  one 
rate  to  another. 

Figure  1  gives  rates  for  LBH,  2PG,  VK,  01  1356  A  and  NI  1493 
A  as  well  as  the  densities  of  0"^(^P)  and  0('s)  and  finally  the  NO 
density  times  lO'®.  The  0*(^P)  and  0(’s)  densities  give  the 
emission  features  2471  A  and  2972  A  when  multiplied  by  their 
appropriate  transition  rates.  The  scaled  NO  density  gives  the 
fluorescence  emission  in  the  y  bands  when  multiplied  by 
appropriately  scaled  g-f actors.  The  results  in  the  figure  refer 
to  high  solar  activity.  Multiple  scattering  calculations  were 
performed  for  01  1356  A.  The  figure  shows  the  emission  rate 
before  and  after  scattering.  The  applied  NO  density  and  needed 
g-factors  come  from  McCoy. Time  dependent  chemistry 
modeling  was  applied  to  VK,  0*(^P)  and  0('s)  . 

Figure  2  shows  volume  emission  rates  for  Rayleigh 
scattering.  The  upper  panel  gives  selected  rates  shortward  of 
about  2500  A  where  pure  absorption  by  O^  and  O^  at  the 
shorter  wavelengths  prevents  the  emission  from  penetrating  into 
the  lower  atmosphere.  The  lower  panel  gives  selected  rates 
longward  of  2500  A  and  shows  the  dramatic  increase  in  emission  at 
low  altitudes  as  O3  absorption  decreases.  The  longest 
wavelength  shown  is  3053  A  since  interest  in  this  report  ends  at 
about  this  wavelength.  It  is  worth  noting,  however,  that  the 
emission  rate  near  the  surface  continues  to  increase  with 
increasing  wavelength  out  to  about  4000  A.  Values  in  excess  of 
2x10^  photons  cm'^  s"'  A''  are  obtained  by  our  model  longward  of 
3200  A. 


7.  SYNTHETIC  SPECTRA 

Separate  spectra  will  follow  for  the  dayglow  and  Rayleigh 
scattering.  Their  combined  emissions  will  be  shown  in  the  next 
section  in  the  form  of  limb  profiles  for  given  wavelength 
intervals  and  FOVs.  Figures  3-7  give  dayglow  spectra  at  high 
solar  activity  for  tangent  altitudes  of  300,  250,  200,  150,  100, 
80,  60,  40,  20  and  0  km.  Each  figure  displays  spectra  for  two 
successive  altitudes.  The  upper  panel  for  any  given  figure 
covers  the  wavelength  range  1200  to  2000  A  while  the  lower  panel 
goes  from  2000  to  3000  A.  The  dominant  features  from  1200  to 
2000  A  are  HI  1216  A,  01  1304  A,  01  1356  A  and  the  system  of  LBH 
bands  starting  at  1273  A  and  continuing  through  the  entire 
wavelength  range.  The  dominant  features  modeled  from  2000  to 
3000  A  are  the  NO  y,  VK  and  2PG  systems  and  the  lines  0'  2471 
and  01  2972.  LBH  also  contributes  at  short  wavelengths  within 
this  region.  We  do  not  show  spectra  at  low  solar  activity 
because  of  the  quantity  of  information  being  presented.  Low  and 


6 


100 


1000  10000  1E  +  5  1E+B 

VOLUME  EMISSION  RATE  (/crr».»3-s) 


1E+7 


Figure  2.  Rayleigh  scattering  volume  emission  rates 
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Figure  3. 
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Dayglow  spectra  for  tangent  altitudes  of  300 
and  250  km  for  high  solar  activity  and  a  solar 
zenith  angle  of  30°. 
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I  high  activity  results  will,  however,  be  presented  in  the  next 

I  section  in  the  form  of  limb  profiles. 

Figure  8  shows  spectra  for  Rayleigh  scattering.  The  chosen 
tangent  altitudes  are  the  same  as  above  except  that  the  first 
four  have  been  deleted.  Above  100  km,  limb  emission  from 
Rayleigh  scattering  is  insignificant  compared  to  the  dayglow. 

I  The  structure  in  the  displayed  spectra  from  2000  to  3500  A  is 

caused  by  O3  absorption  which  peaks  near  2500  A.  The  rapid 
rise  longward  of  2800  A  results  from  the  penetration  of  solar 
photons  into  the  lower  atmosphere  as  discussed  in  the  previous 
section.  This  explains  why  there  is  much  less  rise  at  the  larger 
I  tangent  altitudes  since  for  these,  the  line  of  sight  does  not 

I  penetrate  to  altitudes  where  the  volume  emission  rate  becomes 

large  at  wavelengths  like  3000  A  (see  Figure  2) . 

8 .  LIMB  PROFILES 

The  spectra  shown  above  appear  within  codes  SYNSPEC  and 
RAYSPEC  as  two  dimensional  arrays  versus  wavelength  and  tangent 
altitude.  Limb  profiles  to  follow  were  obtained  by  integrating 
these  arrays  over  specified  intervals  in  both  variables. 
Integration  over  tangent  altitude  (actually  over  look  angle  in 
the  code)  provides  f ield-of-view  averaging  for  either  circular  or 
I  rectangular  FOVs.  Here,  we  chose  circular  FOVs  to  give  the  limb 

resolutions  shown  in  Table  1  (3,  10  and  30  km).  The  respective 
limb  profiles  from  SYNSPEC  and  RAYSPEC  are  placed  on  files  and 
read  into  another  code  where  they  are  joined  for  each  wavelength 
interval.  It  is  these  profiles  which  we  present  here.  We  have 
selected  sixteen  intervals  from  1400  to  3000  A  as  discussed  in 
I  Section  5.  The  next  sixteen  figures  (Figures  9-24)  show  limb 

profiles  for  these  intervals.  The  upper  panel  and  lower  panels 
in  each  figure  respectively  show  the  results  in  megaRayleighs  and 
in  (iWatts  cm’^  sr"'.  Each  panel  contains  six  profiles  for  the 
possible  combinations  of  limb  resolution  and  solar  activity. 

The  differences  from  one  FOV  averaged  profile  to  the  next  are 
not  distinguishable  on  the  scale  chosen  for  displaying  the 
profiles.  Thus,  we  include  Figure  25  which  shows  the  profiles 
from  Figures  15  (2000-2100  A)  and  24  (2900-3000  A)  in  greater 
detail  around  the  100  km  tangent  point  to  illustrate  the  effect 
of  changing  the  FOV  from  3  to  30  km  resolution. 


9.  AZIMUTHAL  EFFECTS 

The  Rayleigh  scattering  results  presented  above  were 
obtained  by  assuming  isotropic  scattering.  Actual  scattering  is 
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Figure  9.  Limb  profiles  over  the  wavelength  interval  1400 
1500  A  for  high  and  low  solar  activity  and  limb 
resolutions  of  3,  10  and  30  km 
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Figure  10.  Similar  to  Figure  9  except  from  1500  to  1600  A 
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Figure  11.  Similar  to  Figure  9  except  from  1600  to  1700  A 
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Figure  12.  Similar  to  Figure  9  except  from  1700  to  1800  A. 

Rayleigh  scattering  begins  to  appear  in  this  band 
but  has  not  been  included  in  this  figure. 
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Figure  13.  Similar  to  Figure  12  except  from  1800  to  1900  A 
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Figure  14.  Similar  to  Figure  12  except  from  1900  to  2000  A 
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Figure  15.  Similar  to  Figure  9  except  from  2000  to  2100  A 
with  the  addition  of  Rayleigh  scattering 
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Figure  17.  Similar  to  Figure  15  except  from  2200  to  2300  A 
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Figure  18.  Similar  to  Figure  15  except  from  2300  to  2400  A 
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Figure  19.  Similar  to  Figure  15  except  from  2400  to  2500  A 
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Figure  20.  Similar  to  Figure  15  except  from  2500  to  2600  A 
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Figure  21.  Similar  to  Figure  15  except  from  2600  to  2700  A 
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Figure  22.  Similar  to  Figure  15  except  from  2700  to  2800  A 
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Figure  23.  Similar  to  Figure  15  except  from  2800  to  2900  A 
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Figure  24.  Similar  to  Figure  15  except  from  2900  to  3000  A 
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described  by  the  Rayleigh  phase  function,  which  in  normalized 
form  is 


where  9  is  the  scattering  angle  and  integration  over  An  gives 
unity.  Figure  26  provides  a  sketch  of  the  geometry  illustrating 
how  0  changes  with  the  sun's  azimuth.  The  largest  deviation  in 
the  calculated  intensity  from  the  isotropic  result  occurs  for 
single  scattering.  In  this  case  the  source  function  5=5^ 
with  the  above  angular  dependence  and  the  angle  0  relates 
exclusively  to  the  sun's  direction  which  introduces  a  dependence 
of  the  sun's  azimuth  angle  for  a  given  elevation  or  solar  zenith 
angle.  The  intensity  for  this  situation  is  easy  to  calculate 
within  our  model  and  will  be  shown  shortly  for  selected 
wavelengths.  As  scattering  increases,  S  becomes  magnified 
relative  to  S^,  approaches  isotropy  for  which  p(0)  =  1/An  and 
thus  ceases  to  have  any  dependence  on  the  sun's  azimuth  angle. 

As  a  first  step  in  examining  the  effect  of  varying  the  sun's 
azimuth  angle,  we  have  separated  S  into  two  parts.  They  are 
and  S-S(j.  As  mentioned  above,  the  Rayleigh  phase  function  is 
applied  to  S^.  Isotropy  is  applied  to  S-S^  which  becomes 
increasingly  accurate  as  the  overall  magnitude  of  S  increases 
relative  to  S^.  The  intensity  equation  then  becomes 

=  linCl  +  cos^ejls^ 

which  has  been  applied  to  obtain  the  results  to  follow. 

Figures  27  -  29  present  limb  profiles  for  Rayleigh 
scattering  based  on  the  above  equation  for  the  wavelength 
intervals  2000-2100  A,  2400-2500  A  and  2900-3000  A.  We  have 
chosen  a  solar  zenith  angle  of  60°  rather  than  30°  here  to 
achieve  a  greater  range  of  azimuth  angles.  The  full  range  is 
30°  to  150°  but  can  be  effectively  taken  from  30°  to  90° 
since  the  phase  function  yields  identical  results  from  90°  to 
150°.  The  top  panels  in  these  figures  show  the  separate 
components  due  respectively  to  and  S-S^.  For  the 
wavelength  intervals  chosen,  we  observe  that  single  scattering 
dominates  over  multiple  scattering.  This  is  true  until  reaching 
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wavelengths  longward  of  3000  A  at  which  point  the  full  atmosphere 
is  available  for  scattering  and  can  produce  S/S^  ratios  in 
excess  of  unity  compared  to  values  between  approximately  . 1  and 
.01  shown  in  Figures  27  -  29. 

The  lower  panels  in  the  figures  show  how  the  upper  panel 
components  combine  for  azimuth  angles  of  30°,  60°  and  90°. 

We  see  that  an  increase  of  nearly  1.75  is  achieved  going  from 
90°  to  30°. 


10.  FUTURE  DIRECTION 

The  following  list  gives  some  of  the  tasks  we  plan  to  pursue 
in  the  near  future.  The  goal  of  these  tasks  is  to  strengthen  and 
further  test  the  model  and  to  make  predictions  available  to  the 
community  for  sensor  development,  systems  architecture  and  scene 
generation. 

1.  Compare  with  satellite  data 

2 .  Add  features 

3.  Investigate  the  need  to  generalize  the  algorithm 

4.  Perform  calculations  both  instrument  specific  and  in 
standard  brightness  units  with  the  latter  being  for 
generation  of  data  bases. 

The  data  obtained  by  AFGL^s  AIRS  experiment  on  Polar  BEAR 
and  by  various  experiments  on  recent  SDIO  missions  will  be 
considered  under  item  1.  Some  of  the  features  we  wish  to  add  are 
01  1304  A,  HI  1216  A  and  various  systems  and  lines  beyond  5000  A. 
The  first  two  of  these  are  presently  in  the  model  for  one  set  of 
specific  conditions  as  noted  in  Section  5.  We  will  investigate 
the  level  of  effort  needed  to  adequately  describe  their  variation 
as  conditions  such  as  solar  activity  and  composition  change.  We 
are  missing  an  important  spectral  region  with  regard  to  Rayleigh 
scattering,  namely  the  Schumann-Runge  band  region  from  about  1750 
to  2000  A.  Approximate  descriptions  of  the  scattering  in  this 
region  will  be  sought.  Guidance  to  generalizing  the  algorithm 
will  come  from  data  analysis  and  comparisons  with  other  models. 

In  particular,  selected  comparisons  are  planned  with  LOWTRAN7  for 
nadir  viewing  to  determine  if  improvements  are  needed  in  the 
Rayleigh  scattering  part  of  the  algorithm  from  UV  to  visible 
wavelengths.  If  so,  these  improvements  will  likely  concern  our 
description  of  the  background  contribution  coming  from  low 
altitudes  and  from  the  surface. 


11.  ACKNOWLEDGMENTS 

This  work  was  supported  by  funding  from  the  Air  Force 
Geophysics  Laboratory  through  Dr.  R.  E.  Huffman  as  part  of  the 
SDIO/SN  ultraviolet  phenomenology  program. 


12.  REFERENCES 

1.  D.  J.  Strickland,  D.  L.  Book,  T.  P.  Coffey  and  J.  A. 

Fedder,  "Transport  equation  techniques  for  the  deposition  of 
auroral  electrons”,  J.  Geophys.  Res.,  2755  (1976). 

2.  D.  J.  Strickland  and  D.  E.  Anderson,  Jr.,  "Radiation 
transport  effects  on  the  01  1356  A  limb  intensity  in  the 
dayglow",  J.  Geophys.  Res.,  88,  9260  (1983). 

3.  D.  J.  Strickland  and  R.  R.  Meier,  "A  photoelectron  model  for 
the  rapid  computation  of  atmospheric  excitation  rates", 

NRL  Memo.  Report  5004,  Naval  Research  Laboratory, 

Washington,  D.C.  2^75  (1982). 

4.  R.  R.  Meier  and  J.-S.  Lee,  "An  analysis  of  the  01  1304  A 
dayglow  using  a  Monte  Carlo  resonant  scattering  model  with 
partial  frequency  redistribution".  Planet.  Space  Sci., 

10,  439  (1982). 

5.  R.  R.  Meier,  D.  E.  Anderson,  Jr.  and  M.  Nicolet,  "Radiation 
field  in  the  troposphere  and  stratosphere:  I  General 
analysis".  Planet.  Space  Sci.,  10,  923  (1982). 

6.  D.  E.  Anderson,  Jr.  and  R.  R.  Meier,  "The  effects  of 
anisotropic  multiple  scattering  on  solar  radiation  in  the 
troposphere  and  stratosphere".  Applied  Optics,  18, 

1955  (1979). 

7.  D.  E.  Anderson,  Jr.,  R.  R.  Meier  and  J.  B.  Kumer,  "Improved 
model  of  Mie  scattering  contribution  to  troposphere  and 
stratosphere  photodissociaton  fluxes".  Applied  Optics, 

19,  1230  (1980). 

8.  D.  E.  Anderson,  Jr.,  "The  Troposphere-Stratosphere  Radiation 
Field  at  Twilight:  A  Spherical  Model",  Planet.  Space  Sci., 

11,  1517  (1983). 

9.  H.  E.  Hinteregger,  "Representations  of  solar  EUV  fluxes  for 
aeronomical  applications",  Adv.  Space  Res.,  1,  39  (1981). 


39 


H.  E.  Hinteregger,  K.  Fukui  and  B.  R.  Gibson, 
"Observational,  reference  and  model  data  on  solar  EUV  from 
measurements  of  AE-E",  Geophys.  Res.  Lett.,  8,  1147 
(1981) . 

A.  E.  Hedin,  "A  revised  thermospheric  model  based  on  mass 
spectrometer  and  incoherent  scatter  data:  MSIS-83",  J. 
Geophys.  Res.,  10170  (1983). 

Atmospheric  Ozone  1985,  World  Meteorological  Organization 
Report  No.  16,  NASA,  Laboratory  for  Atmospheres,  Code  614, 
GSFC,  Greenbelt,  MD  20771  (1985). 

R.  P.  McCoy,  "Thermospheric  odd  nitrogen  2.  Comparison  of 
rocket  observations  with  a  diffusive  transport  chemical 
model",  J.  Geophys.  Res.,  3206  (1983). 

R.  P.  McCoy,  Private  Communication  (1987). 


